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Abstract 


In  this  report  a  systematic  mathematical  method  is  introduced  for  the 
solution  of  problems  involving  two  coupled  modes  in  a  coupled  system  with 
varying  parameters.  These  problems  involve  systems  of  linear  differential 
equations  with  varying  coefficients. 

By  the  use  of  a  linear  transformation  of  the  dependent  variables  and 
a  double  diagonalization  process,  the  coupled  mode  equations  are  reduced  to 
two  decoupled  Riccati  equations.  The  final  form  of  the  general  solution  is 
obtained  in  terms  of  four  varying  coupling  coefficients  and  a  transform 
parameter. 

To  illustrate  some  applications  of  the  method,  solutions  of  two  special 
cases  which  have  been  solved  previously  by  other  workers  are  obt  ined.  The 
solutions  for  a  number  of  special  cases,  in  which  the  varying  coefficients  are 
specified  or  interrelated,  are  also  obtained.  Further  possible  applications 


are  indicated. 


Introduction 


During  the  past  decade,  workers  in  different  branches  of  Physics  and 
Engineering  have  published  a  considerable  amount  of  material  on  the  applica¬ 
tion  of  coupled  mode  theory  to  various  kinds  of  coupled  systems. 1-8  As  far 
as  the  mathematical  theory  is  concerned,  most  researchers  have  considered 
only  special  cases  of  coupled  systems  having,  constant  parameters  or  slowly 
varying  parameters.5*7  In  the  latter  case  approximate  solutions  have  been 
obtained.  Recently,  in  their  study  of  Faraday  rotation  in  a  ferrite  rod,  Huang 
and  Fan  considered  the  case  in  which  two  self-coupling  coefficients  are  equal 
and  one  of  the  mutual  coupling  coefficients  is  equal  to  the  negative-conjugate 
of  the  other.  3  This  report  extends  the  treatment  further  to  more  general  case 
of  varying  parameters,  and  a  systematic  mathematical  method  is  introduced  to 
solve  the  problem  of  two  coupled  modes  in  a  coupled  system. 

In  coupled  mode  theory,  the  behavior  of  a  coupled  system  is  described  in 
terms  of  the  normal  modes  of  the  uncoupled  system.  The  equations  that  charac¬ 
terize  the  behavior  of  two  coupled  modes  takes  the  form: 

da} 

=  A(z)aj  +  B(z)a2  (1-a) 

da  2 

-jg- s  CU)aj  +  D(z)a2  (1-b) 

In  these  expressions,  aj  and  a2  are  mode  amplitudes  of  the  coupled  system,  and 
A(z),  B(z),  C(z),  and  D(z)  are  mode  coupling  coefficients  which  may  be  any 
arbitrary  functions  of  z.  By  the  use  of  a  linear  transformation  of  the  depen¬ 
dent  variables  and  a  double  diagonalization  process,  we  are  able  to  obtain 
general  solutions  of  the  coupled  mode  equations  in  terms  of  the  four  varying 
parameters  A,  B,  C,  D,  and  a  transform  parameter  ho,  which  is  a  particular 
solution  of  a  generalized  Riccati  equation.  Since  the  particular  solutions  of 
the  generalized  Riccati  equation  have  been  widely  studied,  and  numerous  results 
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are  available  in  the  literature  of  non-linear  differential  equations,9-15 
general  solutions  of  equations  (1)  can  be  obtained  for  a  large  number  of 
special  cases.  The  solutions  that  we  have  studied  can  be  applied  to  many 
different  types  of  systems.  What  is  of  interest  here  is  that  the  results  are 
of  such  generality  that  they  allow  the  description  and  classification  of  a 
wide  class  of  devices  in  a  significant  manner. 

Further  applications  of  our  solutions  to  the  structures  of  log-periodic 
antennas,  traveling  wave  antennas,  and  directional  couplers  are  indicated  in 
the  concluding  remarks. 

I.  General  Solutions  of  Coupled-Mode  Equations 

This  section  is  devoted  to  solving  the  two  dimensional  coupled-mode  equa¬ 
tions  by  linear  transformations  of  the  dependent  variables.  In  order  to 
simplify  the  writing,  we  express  the  coupled  mode  equations  in  matrix  form. 


Equations  (1)  in  matrix  form  become 

a'  =  f  a]  a  (3) 

where  the  prime  denotes  differentiation  with  respect  to  the  independent  variable 
Now,  we  introduce  the  following  transformation  of  the  dependent  variables 
in  equation  (3).  11  Let  be  a  continuously  differentiable,  non-singular 

matrix  for  Zj  <  z  <  z2.  Under  the  linear  change  of  variables  a  -*■  g,  where 

l  =  [k]  k  (■:) 


(3)  is  transformed  into 


i’  *  [R]*‘  <[*  3  [R]  -  [R]’>i 


(5) 
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For  our  purpose,  we  choose  [R]  specifically  to  be 


.  r r“  r‘2] .  r eJ 

Lr21  r22  J  L° 


,/D  dz 


It  is  easy  to  show  by  matrix  algebra  that 

i'  =  [s]  i 


where 


CO 


|"sn  si2  "1  _  r 

L s2 1  s22  J  L 


-/ (A-D)dz 


J (A-D)dz 


For  a  second  transformation,  let  us  again  assume  that  CO  be  a  continuously 
differentiable  non-singular  matrix  in  the  same  range  of  z  as  considered  in  the 
first  transformation:  i.e.  Zj  <.  z  <.  z2.  Under  the  linear  change  of  dependent 


variables  g  -*■  y,  where 


i  •  [t] 


(7)  is  transformed  into 


y*  =  [t]-‘  t  C s  3  C T  3  •  [TJ ' >  y 


Let  us  choose  CO  specifically  to  be 


[T]=L  X  ] 


(11) 


Then  (10)  becomes  after  some  mathematical  manipulations 


y’  *  [0? 


(12) 


where 


w  ■  [::: •  [ 


s  12^21  0 


S2lt] 


(13) 


t 12  and  t2i  are  transform  parameters  which  are  particular  solutions  of  the 
following  differential  equations: 


'ir 


r-  -  iMk£4tfgv«v**»  ■<*  ’  • 


[i*-*»  i  jf^iiini 


*12  "  s 1 2  +  s21t?2  =  0 


•21 


S21  +  S12t21 


2  - 


=  0 


(14-a) 

(14-b) 


and  s12»  S21  are  the  matrix  elements  given  specifically  in  (8). 

Since  CO  in  (13)  is  a  diagonal  matrix,  solutions  for  the  y’s  can 
easily  be  found  by  solving  the  decoupled  differential  equations  of  (12). 
They  are 

,  /sl2t21dz 

Yl  =  C1  e 

,  /s21t12dz 
y2  =  c2  e 

where  cj  and  C2  are  constants  of  integration. 

According  to  transformation  (9),  we  have 

,  /s12t21dz 
*1  *  yl  *  t12>,2  *  C1  * 

,  1*12*21** 

*2  *  *21X1  ♦  y2  1  °lt21  * 


(15-a) 

(15-b) 

.  /s2 l*l2dl 

2*12  e 

(16-a) 

,  1*21*12** 

+  c2  e 

(16-b) 

Also,  from  the  transformation  (4),  and  the  expressions  for  the  g's  in  (16), 
we  obtain 

/s12t21dz  .  K,t,,dz 

(17-a) 

—  •  •  a  4. 

/S12t21dZ 


a,  =  c,ru  a 


a2  =  *1*22*21  * 


+  c2r22  e 


/S21,12dz 
/S21 'l2d! 


(17-b) 


Using  the  following  transformation  of  the  dependent  variables  in  (14), 


where  [»]  is  given  as 


rua  “.2-I  r.-/<A-D)dz  0  1 

“  ‘k  J  =  [o  «/<->- J 


(18) 


(19) 
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We  then  obtain  from  (14) 


hj2  -  (A-D)  hj2  ♦  Ch{2  =  B  (20-a) 

h'2l  +  (A-D)  h21  ♦  Bh2j  =  C  (20-b) 


Also,  (17)  becomes  under  the  transformation  of  (19) 


,  /sl2uUhl2dl 

/ s,,u, „h. ,dz 
•  21  12  21 

(21-a) 

1  =  Clrll  *  ♦  C2rllt12 

e 

,  /S 12U22h21dl  , 

/a21Ullh12dl 

(21-b) 

2  "  C lr22t2 1  e  +  C2r22 

e 

Substituting  all  the  elements  of  ,  and  specified  in  (6),  (8), 

and  (19),  we  finally  obtain  the  general  solutions  of  the  coupled  mode  equations 


al  =  C1  6 


a2  *  Clh21  6 


/( A+Bh21)dz  /( D+Chl2)  dz 

Je  *  c2h12  e 

/(A+Bh21)dz  f  /( DtChl2)dz 


+  c2  e 


(22-a) 

(22-b) 


where  h12and  h21  are  transform  parameters  which  are  particular  solutions  of  the 

♦ 

generalized  Riccati  equations  (20-a)  and  (20-b),  respectively.  It  must  be 
noted  that  in  (22)  hJ2  and  h21  should  be  subject  to  the  condition  that 
h 1 2h2 1  "  1  ^  which  follows  from  the  non-singular  requirement  on  transform 
matrix  [t]. 

Realizing  that  parameters  hj2  and  h2i  are  related  through  the  differential 
equations  (20),  we  can  further  eliminate  one  of  these  parameters.  It  can  be 
shown  that  if  hQ  is  a  particular  solution  of  (20-b),  l/hQ  will  be  a  particular 
solution  of  (20-a).  In  order  to  satisfy  the  condition  that  hJ2h21  -  1  i  0, 
l/hQ  and  hQ  cannot  be  used  as  transform  parameters  for  h12and  h2l,  respectively, 
at  the  same  time. 

Let  l/h0  be  a  particular  solution  of  (20-a).  Then  by  Euler's  theorem  con¬ 
cerning  the  generalized  Riccati  equations  (see  appendix  A),  we  find  that 
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h 


0 


exp  /( D-A-2Bh0)dz 

♦  - 

/ B  exp  {/(D-A-2Bh0)dz)  •  dz 


is  also  a  particular  solution  of  (20-b).  Since  the  quantity 

•  exp  /  (D-A-2Bh0)dz 

hj2^21  “  1  -  - 

h0  / B  exp  {/(D-A-2BhQ)dz}  •  dz 

does  not  vanish  in  the  interval  Zj  <.  z  <  z2,  in  general,  l/h0  and  the  above 
expression  are  appropriate  transform  parameters  for  the  system  (22).  The  result 
of  substituting  the  above  particular  solutions  for  the  h's  in  (22-a)  is 


=  c[  exp  /[a  +  B(hQ  + 


exp  /  (D-A-2Bh  )dz 
/ B  exp{/(D-A-2Bh0 )dz)  •  dz 


)J*dz 


♦  c^  exp  /  (D  ~)  dz 


scj  £/  B  exp{  /(D-A-2BhQ)dz)  •  dzj»exp  J  (A  +  BhQ)dz 


+  c2  e*P  /  <D  ♦  iT-)  dz 


From  (20-a),  it  is  found  that 


(23) 


exp  /( D  +  ^-)dz  3  ^  e*P  /(A+ Bhp)dz  (24) 

where  k  is  an  arbitrary  constant.  Then  (23)  is  reduced  to 

*1  s  [C1  /B  exp^  / (D-A-2Bhfl)dz)  •  dz  +  c2J«exp  /(A+Bhfl)dz  (25-a) 

Similarly,  we  have 

a2  ^(cjhg  /  B  exp{/(D-A-2BhQ )dz)  •  dz  +  Cj  exp{/(D-A-2BhQ)dz 

♦  c2hQ  Jexp  /( A+BhQ)dz  (25-b) 

where  hQ  is  a  transform  parameter  which  is  a  particular  solution  of  the  well- 
known  generalized  Riccati  equation  (20-b),  and  ,  c2  are  constants  of 
integration. 
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Following  the  same  procedures,  we  may  obtain  the  solutions  in  a  slightly 
different  form.  They  are 

aj  =  |  Cjh0  +  CjCh^/c  exp{/(A-D-2ChQ)dz}  •  dz 

+  exp{/(A-D-2Ch0)dz)  ]  J • exp/(D+Ch0)dz  (26-a) 

a2  s  ^cl  +  c2  c  exp{/(A-D-2GhQ)dz}  •  dz]  exp/(D+Ch0)dz  (26-b) 

where  h0  is  a  particular  solution  of  the  generalized  Riccati  equation  (20-a), 
and  cj,  C2  are  constants  of  integration. 

Our  general  solutions  of  the  coupled-mode  equations  have  thus  been  obtained 
by  three  successive  linear  transformations  of  dependent  variables  and  a  double 
diagonalization  process.  There  is  a  transform  parameter  involved  in  our  solu¬ 
tions,  which  is  a  particular  solution  of  a  generalized  Riccati  equation.  The 
problem  of  finding  general  solutions  of  the  coupled-mode  equations  is  now  re¬ 
duced  to  that  of  finding  a  particular  solution  of  a  generalized  Riccati  equation. 

The  next  section  will  be  concerned  with  a  discussion  of  the  generalized 
Riccati  differential  equation.  A  method  of  obtaining  particular  solutions 
of  the  generalized  Riccati  equation  is  given  in  Appendix  B. 

II.  The  Generalized  Riccati  Equation 

It  is  customary  to  give  the  name  "generalized  Riccati  equation"  to  any 

equation  of  the  form  . 

^  +  Py  +  Qy2  =  R  (27) 

where  P,  0,  and  R  are  given  functions  of  z. 

This  equation  has  considerable  theoretical  importance,  since  its 
solutions  are  free  from  movable  branch  points  and  can  have  only  movable 
poles.10  It  is  a  special  case  of  the  Abel  equation.  It  is  supposed  that 
neither  R  nor  Q  is  identically  zero.  If  Q  =  0,  the  equation  is  linear; 
if  R  =  0,  the  equation  is  reducible  to  the  linear  form  by  taking  1/y  as 


a  new  variable. 
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It  has  been  shown  by  Euler  that,  if  a  particular  solution  of  the  genera¬ 
lized  Riccati  equation  is  known,  the  general  solution  can  be  obtained  by  two 
quadratures;  if  two  particular  solutions  are  known,  the  general  solution 
is  obtainable  by  a  single  quadrature.  And  it  follows  from  theorems  by 
Weyr  and  Picard  that,  if  three  particular  solutions  are  known,  the  general 
solution  can  be  effected  without  a  quadrature. 

The  equation  (27)  is  easily  reduced  to  a  linear  equation  of  the  second 
order  by  taking  a  new  dependent  variable  u  defined  by  the  equation 

=  1  dlog_u  (28) 

J  Q  dz 


The  equation  then  becomes 


—  at  =  0 

dz2  0  dZ  dz 

Conversely,  if  in  the  general  linear  equation  of  the  second  order 


(29) 


d2u 

dz2 


du 

+  PldJ+P2U  =  0 


(30) 


(where  PQ,  Pj,  P2  are  given  functions  of  z),  one  may  write 

/  Qydz 
u  =  e'  J 


(31) 


From  which  the  equation  defining  y  becomes 


po« 


(32) 


which  is  of  the  same  type  as  (27).  The  complete  equivalence  of  the  generalized 
Riccati  equation  with  linear  equation  of  the  second  order  is  consequently 
established. 

In  the  following  sections,  we  will  apply  our  method  to  a  number  of 
special  cases  for  which  solutions  have  been  previously  established  by  other 


workers . 
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III.  Uniform  Coupling  of  Two  Lossless  Modes  of  Propagation 

Consider  two  waves  with  time  dependence  which  are  weakly  coupled. 

It  can  be  shown  that  the  coupled-mode  equations  for  this  situation  are 

da, 


~dz  -  j6lal  +  C 1 2*2 


(33-a) 


da 


_ 2 

dz 


C21*l  -  ’S2*2 


(33-b) 


The  cJ2  and  c21  are  the  mutual  coupling  coefficients  per  unit  length.  The 
coupling  is  assumed  uniform  over  the  length  of  the  coupler,  so  that  c12  and 
c21  are  independent  of  z.  The  modes  are  assumed  lossless.  For  this  case, 
the  matrix  elements  of  [*]  are  constants.  Let 


Ta(z) 

B(z) 

r-is,  c12-i 

1 — 1 

> 

1 _ 1 

II 

|_C(z) 

D(z)_ 

L  021  -382J 

where  Bj,  B2,  c12,  and  c21  are  independent  of  z. 
From  (34),  (20-b)  becomes 


(34) 


h’  -  j(8l-B2)  h  +  c12h2  =  c21 


A  particular  solution  of  (35)  can  easily  be  found  to  be 


h  -  -  tan 
0  °12  b  cl2 


where 


a  _!l!i 

6d  2 


and 


6b  =  ^  6d2  +  C 1 2 C2 1 


(35) 


(36) 


(37) 


(38) 
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Substituting  (36)  in  (25)  leads  to 

ai  ■  [c. 


c12  "l  -^aZ 

(  -g —  )  sin  8^2  +  c2cos  8bz  e 


( 39-a) 


jBd  8b  1  "j6aZ 

a,  =  I  <ci  +  c2  —  )  sin  8bs  -  <  =  !  -g^  ♦  c2  ^7  >  cos  6bzJ  e 


■[' 


(39-b) 


where 


8a  ' 


8,  ♦  e2 


(40) 


Equations  (39)  can  be  rewritten  in  a  different  form 


■1Y,2 


al  =  A1  e 


+  A2  e 


■]Y25 


(41-a) 


=  1^7  [(WAi 


-DYjZ  -1Y2z-| 

e  +  ^Y2+Bl*  A2  6 


(41-b) 


where 


and 


n  *  6a  +  6b 


(42-a) 


=  6a  -  6b 


(42-b) 


Here  Yj  and  y2  are  normal  mode  propagation  constants  and  Aj,  A2  are  arbitrary 
constants. 

These  solutions  agree  with  those  obtianed  by  other  methods. 


IV.  Coupled-Mode  Description  of  Guided  Wave  Propagation  through  _a  Ferrite  Rod 
A  general  set  of  coupled  mode  equations  for  wave  propagation  in  a  wave¬ 
guide  partially  filled  with  a  gyromagnetic  medium  has  been  derived  by  Huang 
and  Fan.  They  are 


da 
_ e 

dz 

da 

o 

dz 


K  a  -  K  a 
ee  e  eo  o 


K  a  -  K  a 
oe  e  oo  o 


(43-a) 

(43-b) 
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Also 


K  (z)  =  K  (z) 
ee  oo 


(44-a) 


K  (z)  =  -  K  (z) 
eo  oe 


where  the  K's  are  varying  mode  coupling  coefficients. 
In  order  to  solve  this  set  of  equations,  we  write 


[A]  ■ 


A(z)  B(z) 


C(z)  D(z) 


-K  (z)  K  (z) 
ee  eo 


K  (z)  -K  (z) 
eo  ee 


Then  (20-b)  becomes 


h'  -  K  (h2  +  1)  =  0 
eo 


It  is  easy  to  show  that  a  particular  solution  of  (46)  is 


( 44-b) 


(45) 


(46) 


h„  =  tan  /  K  dz 
0  *  eo 


(47) 


Substituting  hQ  in  (25),  we  obtain  the  general  solutions  of  (43): 

al  =  [ci  /  (~Keo)  exp{/2Ke o  tan(/  Keodz)  dz}  *  dz  +  c2  J 

•  exp  /  (Kee  -  Keotan(  J  Kgo  dz)}  •  dz  J 

=  r-  c,/K  sec2(  J  K  dz)dz+  c,"l  •  cos(  /  K  dz)  •  exp  (-/  K  dz) 
I  1  *  eo  1  eo  2  1  '  eo  v  '  ee 

=  |^-  c1tan(/Kgodz)  +  c2  J  •  cos(/Kgodz)  •  exp(-J  K^dz) 

=  f  -  c.  sin(  / K  dz)  +  c,  cos(  /  K  dz)l  •  exp(-  /  K  dz) 

|  1  '  eo  2  *  eo  I  r  *  ee 

dl  -/<Kee+’Keo)dz  d2  -/<Kee-f<eo)dJ 
=  -  e  +  -  e 


(48-a) 
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a.  = 


[-  cx  tan2(  /Kgodz)+  c2  sec2(  /  K^dz) 

+  c«  tan(  /  K  dz)  I  •  cos(  /  K  dz)  •  exp (-/  K  dz) 

*■  4  eo  I  #  eo  '  ee 

=  cos(  J  KeQdz)  +  c2  sin (/  KeQdz)J  •  exp(-/Kgedz) 

d.  -/( k  +  jK  )dz  d„  -/(K  -jK  )dz 
1  '  ee  J  eo  2  '  ee  J  eo 

6 


■<u=- 


n 


] 


where  dj  and  d2  are  arbitrary  constants.  These  expressions  agree  with 
solutions  given  in  Reference  3. 


V.  Case  (i) 


and 


Here  we  consider  the  situation  in  which 

A(z)  =  D(z) 

B(z)  s  C(z). 


For  this  special  case,  the  generalized  Riccati  equation  (20-b)  becomes 


h*  +  B(h2-1)  =  0 


A  particular  solution  of  (50)  is 


h  =  coth  /  Bdz 


The  solutions  which  follow  are  then 


=  ^Cj  /  B  exp{  /  -  2B  coth(  /  Bdz)  dz)  •  dz  +  c2  J 
•  exp  /  {A  +  B  coth(  /  Bdz))  «dz 
=  [Cl  /  Bcsch2(  /  Bdz)dz  +  c2J*(sinh/  Bdz)  •  exp(  /  Adz) 


( 48-b) 

the 


(49-a) 

(49-b) 

(50) 

(51) 


13 


=  (  -  CjCosh  /  Bdz  +  c„sinh  /  Bdz)  exp( J  Adz) 

/  (AtB)dz  /  (A-B)dz 

=  dj  e  +  d2  e  (52-a) 

a2  s  ^cl  cot^2^  /  Bdz)  ♦  Cj  sech2(  /  Bdz)  +  c2  coth(  /  Bdz)} 

•  (sinh  /Bdz)  •  exp(  / Adz) 

/ (A+B)dz  / (A-B)dz 

*  -  d2  e  +  d2  e  (52-b) 

If  B(z)  and  C(z)  are  purely  imaginary,  then  B(z)  =  -  C  (z).  This  situation 
represents  a  passive  coupling  mechanism  in  the  coupled-mode  of  propagation. 


VI.  Case  (ii) 

Case  (ii)  deals  with  the  situation  in  which 


And 


A(z)  =  j  £s(z)  +  Ujn(z)] 
D(z)  *  j  £s(z)  -  y2n(z)J 

B(z)  =  C(z)  =  ju3n(z) 


(53-a) 

(53-b) 

(53-c) 


where  B(z)  and  n(z)  are  any  arbitrary  function  of  z  and  yj,  y2,  y3  are 
arbitrary  constants.  Substitution  of  (53)  in  (20-b)  yields 


h'  +  j(M1+M2)  n(z)  h  +  jy 3  n(z)  h2  =  jy3  n(z)  (54) 

A  particular  solution  of  (54)  is 


hfl  *  j  /it  (^)2  tan(y  /  1  ♦  (j)2  /  n(z)dz)  -  ^ 


(55) 
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The  results  obtained  by  substituting  of  (55)  into  (25)  are 

al  =  £cicos*  /  n(z)dz  +  c2(j  sinA  /  n(z)dz  +  j  cosA  /  o(z)dzf] 


wrw2 


•  exp  j  /  (0(z)  ♦  — r —  n(z)}  •  dz 


( 56-a) 


a2  =  [cl  /  h(z)dz  -  j  cosA  /  n(z)dz) 


M  —  M 

+  c  cosA  /  n(z)dz  ]]  •  exp  j  /  (6(z)  +  -2  6(z)}  •  dz 


( 56-b) 


where 


A 


+  ( 


V^2 

2 


(57) 


When  v3  is  real,  we  have  the  condition  B(z)  =  -  C  (z).  This  is  the  passive 
coupling  case  in  which  the  group  velocities  of  the  waves  are  in  the  same 
direction. 


VII.  Case  (iii) 

Here  we  consider  A,  B,  C,  and  D  to  be  related  in  the  following  manner: 

A(z)  =  j  |^6(z)  +  UjH(z)J  ( 58-a) 

D(z)  =  j  ^B(z)  -  u2n(z)]  ( 58-b) 

C(z)  =  -  B(z)  =  ju3n(z)  ( 58-c) 

For  this  situation,  the  particular  solution  of  (20-b)  can  be  found  to 
be 

hQ  «  -  j  /  (y)2-l  tanh  y3/(y)2-l  /  n(z)dz  +  j  (59) 
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The  general  solutions  for  and  a2  are  then 

al  s  j^cicosh  K  I  n(z)dz  +  c2{jsinh  k  /  n(z)dz  +  ^  cosh  ic  /  n(z)dz) 


_  y  -^2 

•  exp  j(wj+w2)  /  n(z)dz  1  *  exp  j  /  (B(z)  +  —  n(z)}dz  (60-a) 

l2  =  ^Cjijsinh  k  j  n(z)dz  -  j  cosh  x  /  h(z)dz 

+  c2{cosh  x  /  n(z)dz)  •  exp  j(p1+y2)  /  n(z)dzj 


,  uru  2 

exp  j  /  (B(z)  +  — r —  n(z)}dz 


(60-b) 


where 


(61-a) 


w3 


(61-b) 


If  C(z)  and  B(z)  are  purely  imaginary,  C(z)  =  B  (z).  This  is  the  active 
coupling  case  in  which  the  group  velocities  for  the  two  waves  are  in  opposite 
directions. 


VIII.  Case  (iv) 

In  this  case  we  consider  the  relations 

A(z)  *  D(z)  (62-a) 

B(z)  =  t  C(z)  (62-b) 

The  approximation  (62-a)  is  a  fundamental  assumption  commonly  made  when  treat¬ 
ing  a  pair  of  weakly  coupled  modes.  Under  this  assumption,  we  will  neglect 
the  second  term,  which  is  small  compared  with  others,  in  the  generalized 
Riccati  equation. 


If. 


Thus  (20-b)  becomes 


h'  *  B(z)  h2  =  t  B(z) 
h  =  tan  J  Bdz 


(63) 

(64) 


The  expression  (64)  is  a  particular  solution  of  (63)  taking  the  positive  sign, 
and 


h«  =  coth  /  Bdz 


(65) 


is  a  particular  solution  of  (63)  corresponding  to  the  negative  sign.  On 
substitution  of  these  particular  solutions  in  (22),  we  obtain  approximate 
solutions  for  a^  and  a2. 

With  the  upper  sign  in  (62-b),  we  have 

aj  *  j  -c j  £  /  (Bexp  /  (D-A)dz)  •  dz J  cosh  /  Bdz  +  c2sinh  /  Bdz  } 


exp  /  A  dz 


(66-a) 


a2  " 


|  Cjj  £/{Bexp  /  (D-A)dz)  •  dzj  coth 2(  J  Bdz) 

+  exp  /  (A-D)dz  •  sech2  (/  Bdz)  sinh  j  Bd2  +  c2cosh  j  Bdz  j  exp  J  Adz 


(66-b) 


With  the  lower  sign  in  (62-b),  we  have 


ij  »  |  Cj  £/  (Bexp  /  (D-A)dz)  •  dz  J  sin  /  Bdz  t  c2  cos  j  Bdz  ^ 


exp  /  Adz 


(67-a) 


»2  *|^ci{"[/  ^  BexP  /  (D-A)dz)  •  dz  J  tan2  /  Bdz  +  exp  /  (A-D)dz 
•  sec2  /  Bdz  | •  cos  /  Bdz  -  c2sin  /  Bdz  exp  /  Adz 


(67-b) 
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IX.  Concluding  Remarks 

A  systematic  method  has  been  introduced  to  obtain  solutions  of  the  two 
dimensional  coupled-mode  equations  with  varying  coefficients.  Some  special 
cases  for  which  solutions  are  already  known  have  been  given  to  illustrate  the 
application  of  our  method.  The  solutions  of  sections  V-VIII  may  be  applied 
to  physical  problems:  e.g.  directional  couplers,  microwave  antennas,  etc. 

The  coupled-mode  theory  of  propagation  has  been  applied  by  R.  A.  Sigelmann 
and  D.  K.  Reynolds  in  their  study  of  a  traveling  wave  antenna  with  broad  band¬ 
width.8  Because  of  the  systematic  method  of  solving  the  coupled-mode  equations 
developed  here,  it  may  be  possible  to  find  solutions  for  the  log-periodic 
structure  by  the  coupled  mode  approach. 

It  is  also  anticipated  that  the  coupled  mode  theory  as  applied  to  propaga¬ 
tion  problems  will  have  further  application  in  the  design  of  broadband 
directional  couplers  in  which  the  coupling  mechanism  is  of  the  log -periodic  type. 
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Appendix  A 

Euler’s  Theorem  Concerning  the  Generalized  Riccati  Equation: 

It  has  been  shown  by  Euler  that,  if  a  particular  solution  of  the 
generalised  Riccati  equation  is  known,  the  general  solution  can  be  obtained 
by  two  quadratures. 

To  prove  the  result,  let  y0  be  a  particular  solution  of 

^  ♦  Py  +  Qy2  =  R  (A-l) 

and  write  y  *  y0  ♦  The  equation  in  v  is 

^  -  (P  ♦  2Qy0)  v  +  0  =  0  (Pi-2) 

of  which  the  solution  is 

v  exp{  - /(P  ♦  2Qy0)ds}  -  /  Q  exp{-  J(P  +  2QyQ)dz}  •  dz  =  0  (A-3) 

and,  since  b  =  ,  * - *  ,  the  truth  of  the  theorem  is  manifest. 
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Appendix  B 

Particular  Solutions  of  the  Generalized  Riccati  Equation 

Sometimes  one  or  more  particular  solutions  can  be  found  by  inspection 
or  by  a  lucky  guess.  More  formal  procedures  apply  in  dertain  cases. 

I.  The  form  of  the  generalized  Riccati  equation  is 

y *  =  f(x)  +  g(x)y  h(x)y2  (B-l) 

Let  y  =  u(x)/h(x)  and  (B-l)  is  converted  into 

u*  =  F(x)  +  G(x)u  +  u2  (B-2) 

where  F(x)  =  f(x)h(x);  G(x)  =  g(x)  +  .  If  both  F(x)  and  G(x)  are 

polynomials.  See  Sections  (a),  (b),  (c)  in  turn.  If  f(x),  g(x),  and 
h(x)  are  polynomials,  it  will  not  necessarily  be  true  that  G(x)  is  a 
polynomial,  hence  the  methods  of  those  sections  may  not  apply;  see,  how¬ 
ever,  Section  (d)  for  a  possible  procedure. 

(a)  The  equation  becomes 

u'(x)  =  F(x)  +  u2  (B-3) 

and  F(x)  is  a  polynomial  with  G(x)  =  0.  There  are  two  possibilities. 

(i)  The  degree  of  F(x)  is  odd.  There  is  no  polynomial  solution  of 
(B-3);  hence,  none  of  (B-l). 

(ii)  The  degree  of  F(x)  is  even.  Two  possible  polynomial  solutions  may 
exist.  To  find  them,  first  note  that  if  P(x)  were  a  polynomial  of 
even  degree  2n,  then  /  P(x)  could  be  expanded  in  series  of  the  form. 

i  b  i  b. 

.  .  n  n-1  1  2 

P(x)  =  a  x  +  a  ,x  +  ••••  +  an  +  —  +  —  +  •••• 
n  n-1  0  x  x2 

Perform  such  an  expansion  on  /  -  F(x)  but  stop  the  calculation  with 
the  constant  term.  Call  this  result  X(x);  it  is  a  polynomial  of 
degree  n,  if  F(x)  is  of  degree  2n.  The  coefficients  in  the 
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polynomial  could  be  found  by  a  simple  modification  of  the  square- 
root  extraction  method  of  elementary  algebra,  by  the  undetermined 
coefficients,  or  by  expansion  in  a  Maclaurin  series. 

If  the  differential  equation  has  special  polynomial  solutions, 
they  are  given  by 

u  =  t  X(x) 

Test  both  of  them,  for  both,  one,  or  neither  may  satisfy  the  differen¬ 
tial  equation.  If  neither  is  a  solution  of  (B-3),  there  are  no 
polynomial  solutions  of  (B-3)  and  hence  none  of  (B-l). 

(b)  G(x)  i  0.  Let  u  =  w(x)  -  — and  (B-2)  becomes  w2(x)  *  H(x)  +  w 
where  4H(x)  ♦  G2(x)  =  4F(x)  +  2G'(x).  The  equation  in  w(x)  is  similar 
to  (B-3);  hence,  with  slight  modifications  the  procedure  is  very  much 
like  that  in  (a).  Calculate 

Q(x)  =  G2  -  4F  -  2G' 

There  are  two  cases. 

(i)  The  degree  of  Q(x)  is  odd.  There  is  no  polynomial  solution  of 
either  (b-2)  or  (B-l). 

(ii)  The  degree  of  Q(x)  is  even.  Expand  /  Q(x)  as  in  (a)  again,  stop¬ 
ping  with  constant  term  and  call  the  resulting  polynomial  X(x). 
There  are  two  possible  polynomial  solutions  of  (B-2) 

u. .  (£4J« 

Test  both,  for  neither  may  satisfy  the  differential  equation, 
see  (c). 

(c)  Q(x)  is  a  constant.  This  is  the  necessary  and  sufficient  condition 
that  both  solutions  of  (b)  satisfy  the  differential  equation.  There 


are  two  cases. 
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(i)  X(x)  =  k  0.  Introduce  a  new  variable  with  the  relation 

(G-f  k)  ,  1 


u 


2  +  v 


so  that  (B-2)  becomes 


v’(x)  =  1  +  kv 

kx  1 

It  is  separable  with  solution  v  =  ce  +  —  and  the  corresponding 
solution  of  (B-2)  is 


(G+k)  kx  l.-l 

U  =  -  -j-  t  («  ♦  j;-) 


Two  special  polynomial  solutions  of  it  result  with  c  =  0,  ». 

(ii)  X(x)  =  k  =  0.  The  equation  in  v(x)  becomes  v'(x)  +1=0,  which 
is  separable,  with  solution  v  +  x  =  c.  The  solution  of  (B-2)  is 


u  = 


G  1 
2  ~  c-x 


(d)  Polynomial  coefficients: 

Suppose  that  (B-l)  has  the  form 


♦(x)y'  =  f  (x)  +  g(x)y  +  h(x)y2 


(B-4) 


where  all  of  the  coefficients  are  polynomials.  Assume  it  to  have  a 
special  polynomial  solution  y^  =  P.(x)  and  use  a  new  variable  y  =  u  +  yj 
so  that  it  becomes 

4(x)u*  =  F(x)u  +  h(x)u2 


where  F(x)  =  g(x)  +  2h(x)R(x).  If  there  is  a  polynomial  solution  of 
(B-4),  having  the  form  u  =  (x-a)m,  m  >_  1  then  $(x)  will  also  contain 
a  factor  (x-a).  This  property  restricts  the  possible  types  of  poly¬ 
nomial  solutions.  Unlike  the  previous  cases,  where  only  two  possible 
polynomials  may  exist,  the  more  general  solution  can  have  a  larger 
number  of  such  solutions.  See  (e)  for  further  discussion  of  this  case. 
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(e)  Several  polynomial  solutions 

Given  (B-4),  any  specified  number  of  polynomial  solutions  can  be 

constructed.  The  general  solution  can  be  taken  as 

cfj(x)  +  f2(x) 
y  ”  cf3(x)  +  f4(x) 

If  fjf4  -  f2f  3  i  0,  then  (B-l)  can  be  written  in  the  form  of  (B-4) 
with  ♦(x)  =  fjf4  -  f2f3;  f(x)  =  f1f2*  -  fj'f2;  g(x)  =  f2f3'  -  f2f3 

.  V,'  -  f,V;  h(x)  =  f,f„'  - 

Now  suppose  that  the  f^(x)  are  polynomials.  Select  fj(x)  =  0; 
choose  any  polynomials  desired  as  f^x)  and  f4(x).  Assign  n  special 
values  Cj,  c2,  cn  to  the  arbitrary  constant  and  require  that 

f2(x)  =  (c jf 3  +  f4)  (c2f3  +  fi»)  ••••  (cnf3  +  f 4 ) 

The  differential  equation  (B-4)  will  have  (n+1)  polynomial  solutions, 

*2 

y  =  0  and  y  =  -■  y-:—  i  =  1,  2,  •••  ,  n 

l  Cjr3+r4 

The  special  solution  y  =  0  can  be  avoided  and  (n+1)  special  solutions 
retained  if  the  variables  are  transformed  with  the  relation 
u(x)  =  g(x)  +  P(x)  where  P(x)  is  any  polynomial. 

The  considerations  of  this  section  could  be  very  helpful  in  construct¬ 
ing  an  equation  with  a  predetermined  number  of  polynomial  solutions. 

II.  The  following  properties  of  the  generalized  Riccati  equation  sometimes 
apply  in  attempting  to  solve  it. 

(a)  Removal  of  the  linear  term 

The  result,  which  is  similar  to  (B-3)  can  be  achieved  in  three 
different  ways.  It  might  then  be  treated  as  in  I,  but  see  also  (b) 
for  a  special  case. 
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(i)  Let  y  =  u(x)e^X\  $(x)  =  /  g(x)dx  and  (B-l)  becomes 

u'(x)  =  F(x)  +  G(x)u2  (B-5) 

where  F(x)  =  f(x)e  ^X\  g(x)  =  h(x)e^X^. 

(ii)  Let  y  =  u  -  v(x),  v(x)  =  2^X^  .  If  g,  h  are  differentiable  and 
h(x)  t  0,  the  result  is  again  like  (B-5)  but  now  with 

e2 

F(x)  =  f  +  v’  "  =  ' 

(iii)  Let  y  =  u(z)  e^X^;  $(x)  =  /  g(x)dx;  z  =  -  /  he^dx.  In  this  case 
(B-l)  becomes 

u'(z)  =  F(z)  -  u2(z)  (B-6) 

with  F(z)h(z)  =  -  f(z)  e 

(b)  Relation  between  the  coefficients  when  certain  relations  exist  between 
the  coefficient  of  (B-l)  its  solution  may  be  easy  to  obtain. 

(i)  Look  for  two  constants  a,  b  with  |a|  +  |b|  >  0  so  that 

a2f  +  abg  +  b2h  =  0 

If  a  =  0,  then  f(x)  =  0  and  the  equation  is  linear.  If  a  M, 
then  yj  =  b/a  is  a  particular  solution  of  (B-l).  A  simple  case 
arises  iff+g+h=0,  for  then  a  =  b  =  1  and  yj  *  1. 

(ii)  Use  (a) to  remove  the  linear  term  in  (B-l).  Then  if  F(x)  is  propor¬ 
tional  to  G(x)  in  (B-5),  the  result  is  separable.  In  the  original 
variable  and  functions  the  requirement  is  f(x)  =  A2h(x)/exp(2  g(x)dx) 
where  A  is  a  constant  for  proportionality.  The  solution  of  (B-l) 

is  then 

y  =  /  ^  tan (/  /  fh  dx  +  C);  fh  >  0 
if  fh  <  0,  replace  tan  by  tanh  and  insert  a  minus  sign  under  both 


radicals. 
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(iii)  Assume  that  a  special  solution  of  (B-l)  exists  so  that 

2hyj  =  X(x)  -  g(x) 

where  X(x)  is  determined  from  the  relation 
f(x)  =  hyj2  -  X(x)y j  +  y^ 

The  last  equation  imposes  a  severe  restriction  on  the  form 
of  f(x)  but  if  it  is  satisfied,  y^  can  be  found.  A  few 
cases  that  might  be  tested  are 

X  =  0,  4f  =  g2/h  -  2(g/h) '  ; 

X  =  -  h'/h,  4f  =  2(£)"  -  2(g/h) '  -  h(l/h)’2  +  g2/h  ; 

X  =  g  -  2  /"~fh,  2g  =  4  /“fh  ♦  f’/f  -  f'/h  . 

III.  If  all  the  tests  in  I  and  II  fail,  transform  the  equation  into  one  of 
second  order.  Change  the  dependent  variable  by  the  transformation 

yh(x)  u(x)  +  u*(x)  =0;  In  U(x)  +  /  yh(x)dx  =  0 

The  result  i 

u"(x)  +  P(x)u'  +  Q(x)u  =  0 

where  p(x)  s  -  (g  +  h’/h),  Q(x)  =  f(x)h(x).  It  is  linear  and  of  second 
order.  Such  equations  have  been  studied  extensively  and  this  may  be  the 
most  suitable  procedure  unless  the  given  first-order  equation  has  some 
special  property. 

IV.  There  are  fourteen  transforms  listed  below  in  pairs  for  solving  a  certrin 
class  of  Riccati’s  equation.16  The  intermediate  steps  dealing  with 
transforms  and  solutions  of  these  equations  are  avoided.  However,  it 

is  appropriate  to  cite  the  first  case  of  the  transforms  listed  below 
to  explain  the  intermediate  steps  omitted.  Take  transform  1, 


y  =  z/(Pz  ♦  Qz* ). 


(B-7) 
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Substitute  this  into  a  generalized  Riccati's  equation, 

y'  *  a(x)y  +  $(x)y2  ♦  y(x)  (B-8) 

where  a,  $  and  y  are  undetermined  variable  coefficients  and  only  a  can 
become  zero,  typical  of  Riccati's  class  of  equations. 

Substitute  (B-8)  into  (B-7) 

Qz'2  -  P'z2  -  (Q’z'  +  Qz")z  ♦  (Pz  +  Qz')az  +  0z2 

=  (Pz  +  Qz'  ♦  2PQzz' )  (B-9) 

There  are  four  places  where  z2  appears  on  the  other  hand,  there  are  two 
places  where  z'2  appears.  The  key  idea  is  to  eliminate  terms  z2  and 
z'2  in  (B-9),  thus  degenerate  (B-9)  into  a  first-order  linear  differential 
equation.  Because  z'2  appears  in  two  places,  the  unique  choice  for  y(x) 
is  made  first. 


y(x)  = 


QTxT 


(B-10) 


The  coefficients  for  z2  are  collected  to  specify  o(x)  and  8(x), 

p2 

-  P '+  aP  +  0  -  -q*  0  .  (B-ll) 

Note  that  there  are  three  ways  of  specifying  while  retaining  (B-8)  in 
Riccati's  class  of  equations.  If 

a  =  0  (B-12) 


then 


If 


then 


a 


P 

Q 


0'  *  P'  . 


(B-13) 

(B-14) 


(B- 15) 


If 


a  *  (log  P)* 
6 


(B-16) 
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then 


b  - 

S  Q 


( B-17 ) 


Judging  from  t'.iese  developments,  it  is  quite  obvious  that  a  transform 
of  type  . 


y  = 


(B-18) 


Pz  +  Qz' 

does  not  solve  a  Riccati's  equation.  If  (B-18)  is  substituted  in  (B-8), 
there  is  only  one  place  where  z2  appears.  This  forces  P  to  zero  to 
eliminate  the  term  of  z2.  If  P  is  zero  in  (B-18),  y  is  not  related  to 
another  variable  z.  Following  the  above  steps,  the  transforms  listed 
below  are  obtained  and  may  be  applied  for  exact  solutions  for  a  certain 
class  of  Riccati's  non-linear  differential  equations.  Only  fourteen 
transforms  in  pairs  are  listed  below  but  it  is  possible  to  find  many  more 
transforms  in  pairs. 

Transform  1: 

y  = 


solves 


Pz  +  Qz' 


p2 

y'  ♦  (P'  +  —  )y 


1 

Q 


y.  +  h.  .  p.v2  .  I 
y  0  y  Q 

y*  t  (log<P)-y  +  ^  i 


Transform  2: 


solves 


(Pz  +  Qz') 
y  =  - z - 


Y  2  p2 

v '  +  =  p 1  +  £_ 

y  Q  Q 

y *  .  PZ  +  )d  =  p. 

y  Q  Q 


y'  -  (logeP)'y  +  = 


( B-19) 

(B-20) 
( B-21) 
( B— 22 ) 

(B-23) 

(B-24) 

(B-25) 

(B-26) 
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Transform  3: 


solves 


z  t  PQz' 


y'  -  CS(x)  -  ^-]  y  +  6(x)y2  =  ^ 


where  8(x)  is  any  arbitrary  function  of  x. 
Transform  4: 

(z  +  PQz') 

v  =  - 2 - 

3  z 

solves 

y'  +  Cy(x)  “  y  +  pq  =  y(x) 

where  y(x)  is  any  arbitrary  function  of  x. 
Transform  5: 

-  PQz 
y  ”  (Pz  +  Qz' ) 


solves 


*  w2 


..  „  <p  -  Q’>y 


=  p 


Transform  6: 


solves 


y '  +  -q^+  <  Q~ 1 )  *  y 2  =  P 

r  -  (log  Q)’y  +  ^  =  P 

Q2 


(Pz  +  Qz') 
PQz 


y  ♦  py2  =  ^ 

Q2 

y'  -  ^  +  py2  =  (Q'1)1 


y'  +  (log  Q)'y  +  Py2  =  — 
e  Q2 


(B-27) 

(B-28) 

(B-29) 

(B-30  ) 

(B-31) 

(B-32) 

( B— 3  3 ) 
(B-34) 

(B-35) 

(B-36) 

(B-37) 

(B-38) 
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Transform  7: 


y  5 


TpqTTTT 


solves 


Transform  8: 


solves 


y*  +  C(PQ)2  +  <PQ)']y2  =  1 


y’  ♦  PQy  +  (PQ)'y2  =  1 


y'  +  [log  PQ]'y  +  (PQ)2y2  =  1 


-  <PQz  +  z1) 

y  z 


y'  +  y2  =  <PQ)2  +  (PQ)' 


y'  -  PQy  +  y2  =  (PQ)1 


Transform  9: 


solves 


y*  -  [logePQ]'  +  y2  =  (PQ)2 


Pz 

y  =  (Pz  +  Qz') 


y '  -  [8(x)  -  ^]  y  +  8(x)  y2  =  ^ 


where  8(x)  is  any  arbitrary  function  of  x. 
Transform  10: 

(Pz  +  Qz') 


Pz 


solves 


y’  +  Cy(x)  -  y  +  =  y(x) 


(B-39) 

(B-40) 

(B43) 

(B-42) 

(B43) 

(B-44 ) 
(B-45) 
(B-46) 

(B-47) 

( B— 48 ) 


(B-49) 

(B-50) 


where  y(x)  is  any  arbitrary  function  of  x. 
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Transform  11: 


-  Qs 

y  "  (Pz  ♦  Qz') 


solves 


y'  ♦  {[log,  <£)]'  ♦  |  -  2-Mil  |  y  t  B(x)y2 


=  1 


where  0(x)  is  any  arbitrary  function  of  x. 
Transform  12: 


y  = 


_  (Pz  +  Qz') 


Qz 


solves 


V  -(^Cl0ge|r  -?-iliL}yty2  = 


(x) 


where  y(x)  is  any  arbitrary  function  of  x. 
Transform  13: 


PQz 

y  •  <z  ♦  z‘) 


solves 


Transform  14: 


solves 


r  ♦  yJ  5  ^ 

y’  ♦  y  ♦  y2  *  PQ 

y'  -  (logePQ)*  y  ♦  ^  =  PQ 


V  +  PQy2  =  ♦  <^>'3 

y'  -  y  +  PQy2  = 

y'  ♦  (logePQ)'y  +  PQy2  =  ^ 


(B-51) 

(B-52) 


(B-53) 

(B-54) 


(B-55) 

(B-56) 

(B-57) 

(B-58) 

(B-59) 

(B-60) 

(B-61) 

(B-62) 


30 


V.  A  Table  of  Solutions  of  Riccati's  Equations. 


Tabulated  exact  solution  of  y’(x)  +  P(x)y(x)  +  Q(x)y  (x)  =  R(x) 

A  particular  solution 

Required  interrelationship 

(i)  a  WPdx 

R  =  a20  e-j/pdx 

(ii)  1 

/q dx  +  b 

R  =  P 

/Odx  +  b 

(iii)  /Rdx  +  c 

Q  = - L - 

/ Rdx  +  c 

(iv)  -£ 

R  =  -  <§)• 

(V,  f 

0  ■  <£>• 

(Vi)  /~E 

P  =  -  (tn  /"?)’ 

e-/(P-2Q)dx./pe-/(P-2QMx<,x.ii 

(vii)  - 

e-/(P-2Q)dx  +  /pe-/(P-2Q)dx  dx+a 

R  =  0  -  P 

/pe-/(P+20)dx  dxtb_e-/< P+20)dx 

(viii)  - — 

/pe-/(P.20)dx  dx  +bt  e-/(Pt20)dx 

R  =  0  ♦  P 

NOTE:  a,  b,  c  are  arbitrary  constants. 

These  eight  table  entries  represent  a  great  number  of  specialized  cases. 
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